Transverse flow in thin superhydrophobic channels 
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We provide some general theoretical results to guide the optimization of transverse hydrodynamic 
phenomena in superhydrophobic channels. Our focus is on the canonical micro- and nanofluidic 
geometry of a parallel-plate channel with an arbitrary two-component (low-slip and high-slip) coarse 
texture, varying on scales larger than the channel thickness. By analyzing rigorous bounds on the 
permeability, over all possible patterns, we optimize the area fractions, slip lengths, geometry and 
orientation of the surface texture to maximize transverse flow. In the case of two aligned striped 
surfaces, very strong transverse flows are possible. Optimized superhydrophobic surfaces may find 
applications in passive microfluidic mixing and amplification of transverse electrokinetic phenomena. 

PACS numbers: 83.50.Rp, 47.61.-k, 68.08.-p 



Introduction.— Hydrophobic solid surfaces with spe- 
cial textures can exhibit greatly enhanced ("super") 
properties, compared to analogous fiat or slightly dis- 
ordered surfaces [l| . If the liquid follows the topological 
variations of the surface (the Wenzel state), roughness 
can not only significantly increase hydrophobicity, but 
can also lead to a giant drop's adhesion. In contrast, 
when the recessed regions of the texture are filled with 
gas (the Cassie state), roughness can dramatically lower 
the ability of drops to stick and produce remarkable liq- 
uid mobility. At the macroscopic scale, such surfaces 
are "self-cleaning", causing droplets to roll (rather than 
slide) under gravity and rebound (rather than spread) 
upon impact. At the microscopic scale, superhydropho- 
bic surfaces could revolutionize microfluidic lab-on-a-chip 
systems [2|,[3|, which are becoming widely used in biotech- 
nology. The large effective slip of superhydrophobic sur- 
faces [J, [5[ compared to simple, smooth channels [6|-|8[ 
can greatly lower the viscous drag of thin microchan- 
nels and reduce the tendency for clogging or adhesion 
of suspended analytes. Superhydrophobic surfaces can 
also amplify electrokinetic pumping or energy conversion 
in microfluidic devices, if diffuse charge in the liquid ex- 
tends over the gas regions [9rlll|. 
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Superhydrophobic surfaces in nature (e.g. the lotus 
leaf) are typically isotropic, but microfabrication has 
opened the possibility of highly anisotropic textures [l2J . 
The effective hydrodynamic slip [2|, |]J3, LL4| (or electro- 
osmotic mobility [lOf) of anisotropic textured surfaces is 
generally tensorial, due to secondary flows transverse to 
the direction of the applied pressure gradient (or electric 
field [151). I n the case °f grooved no-slip surfaces (Wen- 
zel state) , transverse viscous flows have been analyzed for 
small height variations [16( and thick channels [l7(, and 
herringbone patterns have been designed to achieve pas- 



FIG. 1: Sketch of a thin channel of thickness h(x,y) with 
notation for directions along the plates. In the Wenzel state 
the liquid would conform the solid surface, but in the Cassie 
state it remains free-standing at the top of the roughness 



sive chaotic mixing during pressure-driven flow through 
a microchannel [18l42fJ |. Convection is often required to 
mix large molecules, reagents, or cells in lab-on-a-chip 
devices, and passive mixing by textured surfaces can be 
simpler and more robust than mechanical or electrical 
actuation. In principle, these effects may be amplified 
by hydrodynamic slip (Cassie state) and large amplitude 
roughness (Wenzel state), but we are not aware of any 
prior work. 

In this Letter, we present some general theoretical re- 
sults to guide the optimization of transverse hydrody- 
namic phenomena in a thin superhydrophobic channel. 
We consider an arbitrary coarse texture, varying on scales 
larger than the channel thickness, and optimize its orien- 
tation and geometry to maximize pressure-driven trans- 
verse flow. Our consideration is based on the theory of 
heterogeneous porous materials [21], which allows us to 
derive bounds on transverse flow over all possible pat- 



terns [5J]. 

General considerations.— We consider the pressure- 
driven flow of a viscous fluid between two textured paral- 
lel plates ("+" and "-") separated by h(x,y), as sketched 
in Fig[T] Channel thickness, h, is assumed to vary slowly 
in directions x and y along the plates. We assume a very 
general situation, where sectors of different h are char- 
acterized by spatially varying, piecewise constant, slip 
lengths b + (x,y) and b~(x,y). 

To evaluate the transverse flow, we calculate the veloc- 
ity profile and integrate it across the channel to obtain 
the depth-averaged velocity u in terms of the pressure 
gradient Vp along the plates. As usual for the Hele-Shaw 
cell, the result may be written as a Darcy law 



u= — Vp, 

where the local permeability is: 
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with /3+ = b+/h and 0~ = 
the heterogeneities in h, j3 + 
obtain: 
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(see [5[ for details), we 
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where U,P denote the averages of u,p, respectively. To 
simplify the notation, let G = — (l//i)VP and G = |G|. 
A general inhomogeneous medium is characterized by 
an effective permeability tensor k* with eigenvalues ku 
along x and k± along y, where feii > fcj_ > 0. The vector 
G is applied at an angle a to x. Due to inhomogeneity, 
the velocity U is generally at an angle 9 (0 < 6 < n/2) 
with respect to G. Since ku > fcj_, it is expected that 
U will be preferentially in the direction of x. Letting X 
be the direction of G, and Y the perpendicular direction 
along the plates, we obtain: 



U x = fcii G cos a 



Uu 




FIG. 2: The ratio of velocity components, \Uy/Ux\, as a 
function of A = <Jkn~jkx. 



Therefore, we have transformed our task to optimization 
of A. To maximize Uy /Ux, X should be as large as pos- 
sible (Fig. [2]), i.e. k\\ should be as large, and k± should 
be as small as possible. 

Two-component medium.— In order to proceed fur- 
ther, the analysis is now restricted to a two-component 
anisotropic medium with permeabilities k\,k 2 . Consider 
without loss of generality that < k\ < k 2 < oo. 
The largest possible ku corresponds to the upper Wiener 
bound 
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where <f>\ and 4>2 are the area fractions of the two phases 
with <f>\ + <p2 = 1- The smallest possible k±_ corresponds 
to the lower Wiener bound: 
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Ux = U x cos a + U y sin a = G (k\\ cos 2 a + k± sin 2 a 

) sin a cos a 



A texture satisfying simultaneously both conditions ex- 
ists: it is a configuration of stripes. 
We then have for this texture 
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Our aim is to optimize the texture and the angle a, so 
that the angle between U and G is maximum providing 
the best transverse flow. That is, \Uy/Ux\ should be 
maximum (note that Uy < here). Since at a = and 
7r/2, Uy = 0, these cases are readily eliminated (The 
peculiar case of small \ir/2 — a\ will be treated separately 
below). It is easy to show that 



The surface fraction <pi corresponding to a maximum of 
Uy /Ux can be found from the equation 
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which leads to dX 2 /d<pi = 0, which is satisfied at </>i = 
1/2. This extremum corresponds to a maximum and then 
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A 2 = l 
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with A 2 = feii/A;_L > 1 and t = tana, is at a maximum if 
t = X. The value of the maximum is 
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Defining: 
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the maximum occurs at: 
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and its value is: 
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which increases monotonously with the anisotropy, 
k 2 /ki. It is interesting to note that the preceding analy- 
sis applies to any incompressible, gradient-driven "flow" 
in a two-component medium, not only fluid flow, but also 
electrical conduction (in which case, we have maximized 
the transverse current). 

Superhydrophobic channels.— We now apply these 
results to transverse viscous flow in textured, slipping 
microchannels. We focus first on rough hydrophobic sur- 
faces in the Cassie state, where trapped air bubbles [22|— 
24| can lead to dramatic local slip enhancement. To 
model this, we assume the liquid surface is approximately 
flat, so the local channel thickness h = H is fixed. This 
common assumption (2i 



25| also corresponds to the min- 



imum dissipation [J, |26j. The liquid contacts the solid 
only over an area fraction <p\ of the surface with slip 
length b\ , while the remaining area fraction <p 2 is a free- 
standing gas-liquid interface. As a simple estimate, lu- 
bricating gas sectors of height 8 with viscosity /i g much 
smaller than that of the liquid /x [27[ have a local slip 
length b 2 « 5(/i/iig) ra 50 6, which can reach tens of fim. 
Hydrodynamic slip can also occur at solid hydrophobic 
sectors [28l430l ] , but with b± less than tens of nm [6j-|8j . 

We now consider two cases [5j: (I) one slip wall (/3 + = 
/3; j3~ =0), and (II) equal slip on opposite surfaces (/? + = 
/3~ = /?). Case (I) is relevant for various setups where the 
alignment of opposite textures is inconvenient or difficult. 
Case (II) is normally used to minimize the drag [5|. In 
each case, we have a two-component medium where /3 is 
either /3i or (3 2 . The permeability can now be expressed 
in term of the gap and slip lengths, Eq. ([2]). Then, for 
3 = 1,2: 




case (I) 
case (II) 



(7) 



Since h is constant, the largest k 2 jk\ obviously corre- 
sponds to a largest physically possible fi 2 with smallest 
possible /3i, that is /3i = 0. 

In case (I), approximating the largest possible (3 2 by 
02 — > oo gives k 2 jk\ = 4. We obtain A = k/(k) = 
5/4. The direction of G is a = arctan(5/4) = 51.34°. 
Then \U Y /U X \ -> 9/40 = 0.225. The direction of U is 
a+e = arctan(5/4)-arctan(9/40) = 38.66° (sec Fig. Hi), 
corresponding to a maximum deflection of almost 13°. 




FIG. 3: Results of transverse flow optimization for flat mi- 
crochannels with one superhydrophobic and one no-slip wall 
(case I, left) and two aligned superhydrophobic walls (case 
II, right). In both cases, the thickness of stripes is taken as 
random, to emphasize that the design does not have to be pe- 
riodic. The mean surface fraction of each component is 1/2. 
As an example, the distribution of thicknesses is Gaussian 
with a standard deviation of 0.2. For case II, fa = 100, and 
U x — 300G has been shortened for visibility. 



In case (II), the deflection can be more dramatic, but 
the analysis is more subtle. For f3 2 ^> 1, fe — 6/32, the 
angle a is close to 7r/2. Depending on which of the two 
limits (3 2 — > oo and a —¥ 7r/2 is taken first, the results 
are different. The resolution of this singular perturbation 
problem is to find the significant degeneracy 31] , that is 
the most general limit from which all other cases may be 
obtained. It can be proved here that the significant de- 
generacy is obtained for our optimum. We then calculate 
the following first order approximation: 
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and G is then 9 ~ 7r/2 — 
Note that the flow in this direction 



close to a = tt/2 is large, but is yet O(y r 2 ) smaller than 
the flow that would exist in the X direction for a = 0. 

For completeness, we also apply our results to the Wen- 
zel state, where the liquid is assumed to follow all the 
topological variations of the material. This leads to a 
variable thickness for the liquid domain, but fixed hydro- 
dynamic boundary condition. Therefore, the slip length 
/3 is fixed (possibly zero), but channel thickness h may 
take two values, Hi and H 2 {H 2 > Hi). It is now con- 
venient to define Hi = H on the top of asperities, as we 
defined for the Cassie case. It is easy then to show that 
k 2 /ki = (H 2 /H) 2 . Replacing this value in © shows 
that H 2 /H should be as large as possible. This limit 
is very different from the small surface height modula- 
tions and lubrication geometries considered by Stroock 
et al [18| and suggests that further improvements in pas- 
sive chaotic mixers may be possible with deeper grooves 
and thinner channels. 



Concluding remarks.— A striking conclusion from 
our analysis is that the surface textures which optimize 
transverse flow can significantly differ from those opti- 
mizing effective (forward) slip. It is well known that the 
effective slip of a superhydrophobic surface is maximized 
by reducing the solid- liquid area fraction </>i [4j, |5[ , until 
the Cassie state becomes metastable [12|. In contrast, 
we have shown that transverse flow in thin channels is 
maximized by stripes with a rather large solid fraction, 
4>i = 1/2, where the Cassie state is typically stable. In 
this situation, the effective slip j3* is relatively small [5(, 
and yet the flow deflection is very strong (nearly 7r/2). 

These results may guide the design of superhydropho- 
bic surfaces for robust transverse flows in microfluidic de- 
vices. Applications may include flow detection, droplet 
or particle sorting, or passive mixing. The latter results 
from interactions with side walls, which produce trans- 
verse vortices (due to pressure-driven backflow) and over- 
all helical streamlines, which can be made chaotic for ef- 
ficient mixing by modulating the surface texture in the 
axial direction, e.g. with herringbone patterns 16, 18 
|20[ ■ Compared to the grooved no-slip surfaces with small 
height variations used in prior work, we have shown that 
slipping (Cassie) and highly rough (Wenzel) surfaces can 
exhibit much stronger flow deflection in thin channels, 
which could lead to more efficient mixing upon spatial 
modulation of the texture. 

Another fruitful direction could be to consider trans- 
verse electrokinetic phenomena [15| . e.g. for flow sensors 
or electro-osmotic pumps 



321 ] . It was recently shown 



that flat superhydrophobic surfaces can exhibit tenso- 
rial electro-osmotic mobility |10j : Anisotropy is maxi- 
mized if the Debye screening length Xrj is comparable 
to the texture scale, and the gas-liquid interface is un- 
charged (which, again, does not maximize forward flow); 
the electro-osmotic mobility scales as, M e oc (I + h/Xn) 
in the limit of thick double layers and even thicker chan- 
nels (I < Ad < H) [10J. If a similar relation holds for 
thin channels (H <C L <C AnL then our results for the 
effective b (from k in case II [5J]) suggest that transverse 
electrokinetic phenomena could be greatly amplified by 
using striped superhydrophobic surfaces. 

This research was partly supported by the DFG under 
the Priority programme "Micro and nanofluidics" (grant 
Vi 243/1-3). 
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